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1.  Introduction 

Aluminum  alloys,  as  a  result  of  their  low  mass  density,  high 
strength,  and  high  ductility,  are  of  technological  interest  for  vehic¬ 
ular  applications,  including  aircraft  and  land  vehicles.  The  class  of 
materials  of  present  study  consists  of  bulk  A1  grains  (face  centered 
cubic  or  FCC  crystal  structure)  and  other  phases  that  emerge  as  a 
result  of  incorporation  of  secondary  elements,  either  within  grains 
or  near  grain  boundaries.  The  material  of  interest  here,  A1  2139, 
consists  primarily  of  elemental  A1  and  contains  smaller  amounts 
of  Cu,  Mg,  Ag,  Mn,  and  trace  amounts  of  Si  and  Fe  [3],  The  micro¬ 
structure  may  influence  mechanical  properties  in  a  number  of 
ways.  For  example,  orientations  of  individual  grains  affect  elastic 
compliance  and  yield  strength;  second  phases  influence  disloca¬ 
tion  mechanisms,  and  hence  strength,  ductility,  and  energy  of  cold 
working;  and  precipitates  segregated  at  grain  boundaries  affect 
fracture  strength.  In  Al-Cu-Mg-Ag  alloys  one  important  secondary 
phase  is  the  Q  phase,  with  chemical  composition  Al2Cu  and  ortho¬ 
rhombic  crystal  structure,  which  can  nucleate  as  thin  plates  on 
{111}  planes  and  in  fine  distributions  leads  to  creep  resistance 
and  superior  strength  [14].  A  number  of  other  microscopic  hetero¬ 
geneities  are  also  possible,  including  the  0  phase  (tetragonal  crystal 
structure,  also  of  composition  Al2Cu);  dispersoids  of  Mn;  and  at  the 
single-atom  scale,  point  defects  such  as  interstitials  and  vacancies 
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associated  with  impurities.  Polycrystalline  A1  2139  can  thus  be  la¬ 
beled  a  metal-matrix  composite,  with  the  matrix  phase  itself  a 
composite  mixture  of  anisotropic  grains,  and  with  possible  second¬ 
ary  phases  consisting  of  £1, 0,  and/or  Mn. 

An  increased  understanding  of  mechanisms  affecting  dynamic 
mechanical  performance  of  aluminum  alloys  is  sought.  Relation¬ 
ships  between  composition,  processing  routes,  microstructure, 
strength,  and  ductility  are  needed.  One  goal  of  the  present  work 
is  to  suggest  promising  paths  to  tailor  the  material  at  the  scale  of 
individual  crystals  and  lattice  defect  distributions  in  order  to  im¬ 
prove  performance. 

In  the  present  paper,  a  new  constitutive  model  capturing  defor¬ 
mation  of  individual  anisotropic  crystals  with  effects  of  precipitates 
and  impurities  is  developed,  as  presented  in  Section  2.  A  unique  fea¬ 
ture  of  the  model  not  included  in  previous  work  on  metals  [4-6]  is 
explicit  consideration  of  the  residual  elastic  volume  change  result¬ 
ing  from  the  energy  of  local  stress  fields  of  lattice  defects  [28,12,35], 
More  precise  values  of  thermoelastic  properties  not  considered  in 
Clayton  [6]  are  also  used  here,  following  Thomas  [32].  The  model 
is  implemented  in  a  dynamic  finite  element  code;  a  representative 
simulation  of  grain  morphology  presented  in  Section  3  indicates 
heterogeneous  localized  deformation  resulting  from  spatial  varia¬ 
tions  in  properties  within  the  microstructure.  In  Section  4,  a  new 
macroscopic  model  accounting  for  homogenized  behavior  of  the 
polycrystalline  metal-matrix  composite  is  presented,  capturing 
strain-dependent  variations  in  stored  energy  of  defects  in  the 
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balance  of  energy  and  residual  volume  changes  in  the  equation  of 
state.  The  coarse-grained  model’s  parameters  are  obtained  via  cal¬ 
ibration  to  results  of  the  crystal  plasticity  model,  which  itself  is  too 
computationally  intensive  to  enable  large  scale  ballistic  simula¬ 
tions.  Instead,  macroscopic  finite  element-smooth  particle  dynamic 
simulations  incorporating  the  homogenized  material  behavior  are 
used  to  provide  quantitative  assessments  of  possible  benefits  of 
texturing  and  amplification  of  energy  storage  mechanisms  on  bal¬ 
listic  performance,  as  discussed  in  Section  5.  Such  information 
may  be  useful  to  alloy  developers  seeking  processing  paths  for 
engineered  materials  at  nano-  and  microscales  that  offer  optimal 
performance  in  structural  applications. 


where  on  the  left,  p  =  JEp  is  the  mass  density  of  the  elastically  un¬ 
loaded  material,  p  is  the  mass  density  of  the  deformed  material,  and 
JE  =  detFE.  Components  of  elastic  strain  are  2 EEap  =  FE^  5abFEb  -  6ap, 
with  Kronecker’s  delta  denoted  by  <5,  and  summation  is  invoked 
over  repeated  Greek  and  lower-case  Roman  indices.  Second-  and 
third-order  elastic  coefficients,  C^xS  and  may  generally  de¬ 

pend  on  temperature,  as  may  the  effective  shear  modulus  p.  Alumi¬ 
num  single  crystals  exhibit  three  independent  second-order  elastic 
coefficients  and  six  independent  third-order  coefficients  [32].  Also 
in  (5),  k  and  Q  are  dimensionless  scalars,  c  is  the  specific  heat  per 
unit  mass,  and  0o  is  a  scalar  temperature  at  which  the  thermal  en¬ 
ergy  vanishes.  The  scalar  internal  state  variable  c  is  defined  by 


2.  Model  for  anisotropic  crystal  mechanics 


£  =  by/p^, 


(6) 


The  present  treatment  of  geometrically  nonlinear,  single  crystal 
elastoplasticity  extends  previous  theory  [4,5]  to  account  for  free 
energy  accumulation  associated  with  interactions  between  dislo¬ 
cations  and  precipitates  or  interstitials  as  well  as  residual  elastic 
volume  changes.  Standard  notations  of  modern  continuum 
mechanics  are  used.  Boldface  type  is  used  for  vectors  and  tensors, 
and  italic  font  is  used  for  scalars.  Juxtaposition  of  boldface  quanti¬ 
ties  implies  matrix  multiplication.  The  summation  convention  ap¬ 
plies  over  repeated  indices.  For  simplicity,  Cartesian  coordinates 
are  used  throughout. 

The  local  deformation  gradient  F  satisfies 

F  =  dx<p  =  FEF('j1/3Fp.  (1) 

where  X  are  material  coordinates,  x  =  cp(X,  t)  are  spatial  coordi¬ 
nates  depending  on  time  t.  FE  accounts  for  mechanical  elastic  defor¬ 
mation  and  rigid  body  rotation,  F*  accounts  for  thermal  expansion 
or  contraction,  and  Fp  accounts  for  dislocation  glide.  The  scalar}  ac¬ 
counts  for  residual  elastic  volume  changes  associated  with  nonlin¬ 
ear  elastic  behavior  and  eigenstresses  of  defects  [28,12].  The 
spatial  velocity  gradient  following  from  (1)  is 

L  =  dx\  =  FEFE  1  +  aTn  +  (1/3  )J]-H  +  FEF"LPF''-1FE-1,  (2) 

where  the  superposed  dot  and  -1  denote  a  material  time  derivative 
and  matrix  inversion.  For  cubic  lattices,  thermal  deformation  is  iso¬ 
tropic.  The  coefficient  of  thermal  expansion  aT  may  depend  on  tem¬ 
perature  0,  and  1  is  the  unit  tensor.  Deviatoric  plastic  deformation 
satisfies  the  usual  rate  kinematics 

Lp  =  fpFp-’  =  y(“)s'a)  ®  m<a),  (3) 

a=l 

where  s^1  and  denote  the  reference  slip  director  and  slip  plane 
normal  for  glide  system  a,  with  n  the  total  number  of  systems,  and 
■yw  is  the  shearing  rate  on  system  a.  Slip  directors  and  plane  nor¬ 
mals  are  mapped  from  the  reference  to  deformed  configuration 
via  the  thermoelastic  deformation  as  follows: 


s(“)  =  FEF"j1/3s 


(*) 

0  ’ 


m(“>  =  m‘)°°J-1/3F"-1FE-1. 


(4) 


In  FCC  crystals,  usual  Burgers  vectors  for  mobile  dislocations  are  of  the 
type  a/2(110),  with  magnitude  b=  \/2a,  where  a  is  the  lattice 
parameter  of  the  conventional  unit  cell.  Glide  planes  are  of  orientation 
{111},  leading  to  12  independent  slip  systems  of  type  (1 1  0){1 11}. 

The  Helmholtz  free  energy  per  unit  mass  is  denoted  by  \jt.  The 
specific  free  energy  density  of  the  crystal  is  written 


P t 


—  Fe  (PaPXs  FE  -\-  —  Fe  (T0CPX5e<I>Fe  fe 
2  aft  Xs  “r  g  a/}^  Xs  £<f> 

second-order  elasticity  third-order  elasticity 


-pcB\og(B/80), 

V  ✓ 


(5) 


where  pT  is  the  total  line  length  per  unit  volume  of  dislocations,  in 
general  both  mobile  and  immobile,  associated  with  stored  energy  of 
cold  working  [26]  and  strain  hardening  during  plastic  flow.  Thermo¬ 
dynamic  arguments  [4,5]  lead  to  the  stress-strain  relations 

^  =JEFEa'«aabFEb-'e  =  p^r-  =  C^%  +  lc«W£%sE%,  (7) 

with  S  the  elastic  second  Piola-Kirchhoff  stress  and  a  the  Cauchy 
stress.  The  local  energy  balance  in  the  spatial  configuration  is 


pcB 

temperature  change 


PIT1 


a=l 


residual  volume  changes 


dissipation  from  slip 

_nf  ^  P  9  ^ 

f  l  dz  dBdZ 


energy  of  lattice  defects 


+  8 


pd  i/i  It  dp 
Odd  Ee  <9Ee 


thermoelastic  coupling 


:  EE  -  <9x-q  , 

_ ✓  heat  conduction 


(8) 


where  tw  =  a  :  s(a)  ®  m|a)  is  the  resolved  shear  stress  on  system  a,  p 
is  the  Cauchy  pressure,  and  q  is  the  heat  flux  satisfying 

q  =  ~kdx8,  (9) 


where  k  is  the  scalar  thermal  conductivity.  Symbols  ■  and  :  in  (8)  de¬ 
note  contraction  over  one  and  two  sets  of  indices,  respectively. 

Following  the  generic  Arrhenius  form  for  thermally  activated 
kinetics, 

yM  =  y0exp[-Cw/(M)],  (10) 

where  y0  is  a  constant  with  dimensions  of  1  /£,  kB  is  Boltzmann’s 
constant,  and  Gw  is  the  Gibbs  activation  energy  that  depends  upon 
the  state  variables  [4,5,33]: 

Cw  =  -mkBe[log(TW/goa))  ~  rlog(0/0b)].  (11) 

Here  m  and  r  are  dimensionless  parameters  denoting  rate  and  tem¬ 
perature  sensitivity,  respectively,  and  t(a>  =  JE t(oc).  The  larger  the  va¬ 
lue  of  in,  the  less  rate  sensitive  the  flow  stress.  A  negative  value  of  r 
leads  to  thermal  softening  from  increased  dislocation  mobility  with 
increasing  temperature.  The  transient  slip  resistance  gf/1  = 
g}a)  +  Ag,a)  arises  from  the  sum  of  (i)  the  initial  distribution  of  de¬ 
fects,  lattice  friction,  and  viscous,  phonon,  and  electron  drag  via 
constant  gf1 ;  and  (ii)  evolving  barriers  to  slip  such  as  accumulated 
forest  dislocations  via  state  variable-dependent  contribution 
A g{<x)(Z,  t).  The  evolution  of  slip  resistance  at  reference  temperature 
0o  is  dictated  by  a  hardening-minus-dynamic-recovery  relation: 


go”0  =  A  J2  q;\jw |  -  Bgtf1  1 7(W I ■ 

p=t  0=1 


^  =  <5“  +  g(l-<5“), 


(12) 


lattice  defects 


thermal  energy 
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where  q  is  the  latent  hardening  ratio  equal  to  unity  in  what  follows 
[22],  such  that  all  n  systems  harden  equally.  In  (12),  A  and  B  are 
constants.  The  constant  g}a>  is  independent  of  strain  rate  and 
temperature  and  reflects  the  initial  yield  stress  at  0  =  0o.  Thus, 
=  d(Ag(a>)/df,  where  gj,a)  is  evaluated  according  to  (12).  Possible 
influences  of  rate  and  temperature  on  the  evolution  of  hardness  g(0x> 
are  apparent  upon  substitution  of  the  slip  rates  from  (10)  and  (11) 
into  (12).  The  zero  subscript  on  the  notation  gj,”1  is  not  necessary, 
but  it  is  used  for  consistency  with  previous  work  [4,5,33].  The  origin 
of  this  subscript  follows  from  the  argument  that  at  0  =  0o,  temper¬ 
ature  does  not  affect  the  apparent  yield  stress  entering  the  flow  rule 
(10)  resulting  from  the  form  of  the  activation  energy  in  (11). 

The  following  relationship  between  hardening  and  the  accumu¬ 
lated  dislocation  line  density  is  used  [29]: 

A g(”>  =  dtpb^.  (13) 

Scalar  proportionality  factor  a  accounts  for  dislocation  interactions. 
The  residual  elastic  energy  per  unit  volume  from  dislocations  Es  is 
estimated  as 


Es  =  ^  KQpt2  =  ^KQpb2pT  = 


pb  R 
-4nl0SRc 


QpT 


(14) 


The  first  two  equalities  in  (14)  offer  a  phenomenological  represen¬ 
tation  of  the  stored  energy  of  the  local  stress  fields  of  dislocations. 
Such  a  linear  relationship  between  the  scalar  dislocation  density 
and  stored  energy  has  been  described  elsewhere  [23,4,5];  in  this 
relationship,  the  product  kQ  is  simply  a  proportionality  constant 
to  be  calibrated  from  experimental  measurements  on  stored  energy 
of  cold  work  [26],  The  term  following  the  final  equality  includes  the 
well-known  logarithmic  energy  per  unit  length  of  a  screw  disloca¬ 
tion  in  an  isotropic  body  of  infinite  extent  [13].  Multiplying  this  en¬ 
ergy  per  unit  length  by  the  product  of  the  line  length  per  unit 
volume  of  dislocations  and  the  scaling  factor  Q  then  gives  the  total 
energy  per  unit  volume.  The  relationship  R/Rc  =  exp(2nK)  is  ob¬ 
tained  by  straightforward  algebra.  For  straight,  non-interacting 
screw  dislocation  lines,  a  typical  estimate  is  kQ  ~  2  [13],  Contribu¬ 
tions  of  k  to  deviations  from  this  approximate  value  (i.e.,  kQ^=2) 
account  for  elastic  anisotropy,  core  effects,  elastic  nonlinearity, 
edge  and  mixed  character  of  dislocations,  and  dislocation  line 
curvature.  The  parameter  Q  is  described  in  more  detail  below.  The 
volume  change  associated  with  the  residual  elastic  energy  of  (14) 
is  [28,12] 

J  =  1  +  /t-1] dp/dp  -  p/K)Es  =  1  +  (KQb2 pT/2)(dp/dp  -  p/K), 

(15) 


Furthermore,  from  (5)  and  (13), 


F  ddd^  &pn  dd  s 
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Pt  : 
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Substituting  (17)  and  (18)  into  (16)  then  yields 


(18) 


P=  i  -KQ^p+e^j 


( dipn )2  ^2  T(a)y(a) 


(19) 


As  is  clear  from  (8)  and  (19),  when  k  is  held  fixed,  values  of  Q  >  1 
amplify  the  stored  energy  associated  with  internal  variable  {  and 
equivalently,  the  residual  elastic  energy  engendered  by  microscopic 
stress  fields  arising  from  the  dislocation  density  pT.  Such  amplifica¬ 
tion  would  manifest  physically  in  Al-Cu-Mg-Ag  alloys  from  the 
presence  of  second-phase  particles  of  0,  or  Mn  that  can  act  as  lo¬ 
cal  stress  risers  [2]  or  barriers  to  dislocation  motion.  The  increase  in 
stored  energy  would  be  attributed  to  both  the  elastic  energy  of  the 
deforming  second  phases  and  misfit  stresses  [17],  as  well  as  the 
elastic  interaction  energies  of  dislocations  that  may  build  up  at 
grain  and  phase  boundaries.  While  K  accounts  for  the  self-energy 
of  dislocations  as  discussed  following  (14),  Q  accounts  for  interac¬ 
tion  energies  of  dislocations  with  second  phases,  point  defects  or 
impurities,  and  other  dislocations.  If  the  dislocation  density  van¬ 
ishes  in  a  reference  configuration  of  minimum  free  energy  (e.g.,  a 
perfect  lattice)  then  the  product  kQ  must  be  positive  to  ensure  that 
the  free  energy  (5)  increases  with  increasing  defect  content.  Values 
Q  <  1  are  permissible  and  would  reduce  the  stored  energy  associ¬ 
ated  with  lattice  defects  when  k  is  held  fixed. 

Parameters  entering  the  single  crystal  model,  many  specific  to 
alloy  A1  2139,  are  compiled  in  Table  1.  Because  in  applications  of 
present  interest,  slip  takes  place  before  large  deviatoric  elastic 
strains  can  arise,  third-order  elastic  constants  in  (5)  are  repre¬ 
sented  by  the  simple  form  [5,7] 


(20) 


where  Kt  >  0  accounts  for  the  increase  in  apparent  bulk  modulus 
with  increasing  pressure.  Elastic  constants  and  their  pressure  and 
temperature  derivatives  are  obtained  from  Thomas  [32].  The 
remaining  parameters  in  Table  1  are  obtained  from  the  plasticity  lit¬ 
erature  or  calibration  to  experimental  data  for  A1  2139  [3],  with  a 
random  texture  as  explained  in  Clayton  [6[.  The  large  value  of  m  im¬ 
plies  rate  insensitivity.  A  range  of  values  of  1  sg  Q  ^  8  is  considered 
to  investigate  effects  of  stored  energy  on  the  thermomechanical  re¬ 
sponse,  as  explained  in  detail  in  Section  4. 


3.  Mesoscale  simulation 


where  for  a  cubic  crystal  with  second-order  constants  Cn,C12,  and 
C44,  the  appropriate  bulk  and  shear  moduli  I<  and  p  satisfy 
3 1<  =  Cn  +  2Ci2  and  2 p2  =  €44(011  -  C12).  The  elastic  coefficients 
and  the  pressure  derivative  of  the  shear  modulus  in  (15)  are  all 
measured  at  the  stress-free  reference  state.  The  Taylor-Quinney 
factor  p  dictating  plastic  work-to-heat  conversion  in  adiabatic  pro¬ 
cesses  [10,31,26]  is  defined  as 


P  =  (  J2  ^¥a)  -  P(d 'I'/dZ  ~  0dV/<90(9{){ 


(16) 


such  that  1  -  p  is  the  ratio  of  the  rate  of  stored  residual  elastic  en¬ 
ergy  to  dissipation  resulting  from  dislocation  glide.  Manipulating 
(6)  followed  by  time  differentiation  gives 


{  = 


pdtp/dj 
KpQ  apn  ^ 


—  V  Ag<,z), 

1  in  ^ 
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2 ctpn^ 
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A  direct  numerical  simulation  of  individual  anisotropic  grains  is 
performed  first  to  provide  insight  into  the  high  rate  mechanical 
response  of  the  alloy  at  the  scale  of  the  microstructure. 


Table  1 

Parameters  for  crystal  model. 


Parameter 

Value 

Parameter 

Value 

Cn  [GPa] 

108 

«r  [1/K] 

2.3(10r5  +  1.7(10rsAfl 

C12  [GPa] 

60.4 

c  |J/kgI<] 

900  +  O.466A0 

C44  [GPa] 

28.3 

gj/1  [GPa] 

0.155 

dCn/dB  [MPa/K] 

-35.1 

q 

1.0 

9C,2/afl  [MPa/K] 

-6.7 

a 

0.3 

dCu/de  [MPa/K] 

-14.5 

b  [nm] 

0.286 

dp/dp 

1.75 

yo  [i/s] 

0.001 

/Cl  [GPa] 

36.7 

m 

50 

Po  [kg/m3] 

2700 

n 

12 

K 

10 

A  [GPa] 

0.48 

Bo  [K] 

300 

B 

1.5 

01  [«] 

934 

r 

-1.1 

k  [W/miq 

237 

Q 

1.0,  4.0,  8.0 
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Representative  results  are  shown  in  Fig.  1 .  Properties  are  chosen  to 
correspond  to  an  A1  2139  alloy  of  the  type  described  in  Table  1, 
with  Q  =  1.  As  will  be  demonstrated  later  in  Section  4,  this  leads 
to  an  average  value  of  the  heat  dissipation  fraction  fi  on  the  order 
of  0.90-0.95  when  considered  over  a  range  of  strains  from  zero  to 
unity.  Such  a  prediction  is  in  general  agreement  with  experimental 
observations  of  Farren  and  Taylor  [10],  Havner  [11],  and  references 
therein,  who  noted  for  commercially  pure  polycrystalline  alumi¬ 
num  an  average  value  of  p  ~  0.93.  In  the  present  model,  effects 
of  inclusions  and  precipitates  are  treated  implicitly  by  the  choice 
of  material  parameters,  meaning  that  the  second  phase  particles 
are  not  modeled  explicitly.  In  the  present  context,  the  impurities 
are  too  small  to  resolve  efficiently  using  continuum  finite  ele¬ 
ments:  Mn  dispersoids  are  of  diameter  0.1  -0.5  pm  [3],  and  C2 
platelets  are  of  thickness  less  than  0.1  pm  [14].  The  size  of  the 
computational  domain  is  1  mm  x  1  mm,  consisting  of  100  grains 
of  an  average  size  of  1 00  pm,  with  a  random  initial  texture.  Polyg¬ 
onal  grain  shapes  are  created  using  a  Voronoi  technique  [9,25] 
with  selective  mesh  refinement  in  the  vicinity  of  grain  boundaries 
[4,5].  The  finite  element  method  with  explicit  dynamics  and  heat 
conduction  is  used  in  the  calculations,  as  described  in  [4,5], 

The  boundary  conditions  are  plane  strain,  with  an  applied 
stretch  in  the  vertical  direction  resulting  from  a  velocity  of  10  m/s 
applied  to  the  upper  surface  and  the  bottom  surface  fixed  in  the 
vertical  direction.  An  initial  velocity  gradient  component  in  the 


direction  of  stretch  of  magnitude  104/s  is  assigned  throughout 
the  domain  to  minimize  inertial  effects.  These  boundary  conditions 
provide  an  average  applied  tensile  strain  rate  of  104/s.  Plane  strain 
boundary  conditions  are  used  for  the  two-dimensional  simulations 
because  they  are  deemed  more  representative  (than  plane  stress, 
for  example)  of  conditions  experienced  by  a  polycrystalline  volume 
element  embedded  within  a  much  larger  sample  of  material  [38,4]. 
Special  force  boundary  conditions  applied  along  the  lateral  edges 
constrain  these  edges  to  remain  vertically  straight  and  parallel, 
yet  free  to  contract  inwardly  due  to  the  Poisson  effect  as  the  mesh 
undergoes  tensile  stretch.  The  lateral  edges  are  free  of  shear  trac¬ 
tion.  A  uniform  acceleration  in  the  direction  of  contraction  is  as¬ 
signed  to  all  nodes  of  a  lateral  edge,  with  the  value  of  this  nodal 
acceleration  computed  for  each  explicit  time  step  by  dividing  the 
total  reaction  force  along  the  edge  by  the  total  mass  of  the  nodes 
comprising  that  edge.  A  similar  methodology  was  used  in  explicit 
simulations  of  polycrystalline  microstructures  by  Zhou  et  al.  [38] 
for  imposing  periodic  shear  deformations  by  averaging  forces  over 
periodic  pairs  of  nodes.  The  boundary  constraints  are  intended  to 
instill  the  mesh  with  a  global  deformation  mode  that  could  be  rep¬ 
resentative  of  a  polycrystalline  volume  element  embedded  within  a 
much  larger  sample.  For  example,  spatial  coordinates  of  the  bound¬ 
aries  move  in  an  affine  manner  according  to  x  =  FX,  where  compo¬ 
nents  F\  and  F\  measure,  respectively,  the  average  lateral 
contraction  and  average  applied  stretch,  F\  =  1  in  plane  strain, 
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Fig.  1.  Simulated  Al  polycrystal  with  random  grain  orientations  at  applied  tensile  strain  of  8%:  (a)  effective  stress,  (b)  effective  plastic  strain,  (c)  temperature,  (d)  dissipation 
fraction,  (e)  dislocation  density,  (f)  residual  volume  change. 
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and  FaA  =  0  for  a^A.  These  affine  displacement  conditions  are  anal¬ 
ogous  to  a  macroscopic  deformation  gradient  F  imposed  on  a  single 
linear  finite  element  containing  a  representative  volume  of  micro¬ 
structure.  Although  the  assumption  of  plane  strain  may  inhibit 
the  formation  of  localized  deformation  modes  such  as  shear  band¬ 
ing  and  necking,  the  plane  strain  condition  has  been  used  exten¬ 
sively  by  a  number  of  authors  to  study  shear  localization 
phenomena  in  the  context  of  crystal  plasticity  or  viscoplasticity 
[21,38,1].  The  constraint  imposed  on  the  lateral  edges  of  the  mesh 
in  the  present  work  may  also  inhibit  shear  localization;  for  exam¬ 
ple,  necking  of  the  sample  would  be  possible  if  the  lateral  edges 
were  completely  free  of  traction.  Thus,  the  combination  of  plane 
strain  and  linear  displacement  conditions  used  here  is  expected 
to  provide  a  stiffer  response  than  other  possible  choices  of  bound¬ 
ary  conditions  such  as  plane  stress  or  applied  traction  conditions 
(e.g.,  null  tractions  along  the  lateral  edges).  The  boundary  condi¬ 
tions  used  here  are  sufficient  to  provide  some  insight  into  the  phys¬ 
ics  of  shear  localization  in  a  representative  polycrystalline  sample 
via  examination  of  the  field  variables  (Fig.  1 ),  but  a  more  thorough 
investigation,  reserved  for  future  work,  would  also  consider  of  other 
kinds  of  boundary  conditions.  Three-dimensional  elastic  and  plastic 
deformations  are  permissible,  with  the  full  number  of  potentially 
active  slip  systems,  so  long  as  total  deformation  remains  planar, 
as  assumed  in  a  number  of  previous  works  [8,4,5,33].  A  uniform  ini¬ 
tial  temperature  of  300  K  is  assigned  along  with  null  heat  flux  con¬ 
ditions  normal  to  the  external  boundary.  Several  mesh  densities 
were  investigated,  and  to  the  strain  levels  shown  in  Fig.  1,  the  re¬ 
sults  are  insensitive  to  the  mesh  size. 

Results  in  Fig.  1  correspond  to  an  8%  applied  strain,  that  is, 
t  =  8  ps.  Shown  in  Fig.  1(a)  is  the  scalar  effective  stress,  found  by 
a  =  ((3/2)<r'  :  <t')1/2,  with  a'  the  deviatoric  part  of  the  Cauchy  stress. 
Most  grains  support  stresses  on  order  of  0.5  GPa,  but  owing  to  orien¬ 
tation  and  elastic-plastic  anisotropy,  some  grains  are  stiffer  or  more 
compliant  than  others.  Shown  in  Fig.  1(b)  is  the  cumulative  effective 
plastic  strain  found  by  sP  =  f(( 2/3)Dp  :  D p)1/2dt,  where  Dp  is  the 
symmetric  part  of  Lp  in  (3).  Several  bands  of  localized  deformation 
are  apparent,  most  notably  that  on  the  lower  right  side  of  the  domain, 
in  which  the  plastic  strain  approaches  0.5.  As  is  clear  from  Fig.  1(c), 
regions  of  intense  strain  correlate  with  regions  of  greatest  tempera¬ 
ture  rise,  with  temperatures  on  the  order  of 450K  attained  in  portions 
of  the  band.  As  shown  in  Fig.  1(d),  the  dissipation  fraction  p  correlates 
positively  with  temperature  rise  and  plastic  strain  localization.  Note 
the  inverse  color  scale  used  in  Fig.  1(d),  such  that  red  areas  corre¬ 
spond  to  an  increased  rate  of  energy  storage.  This  correlation  is  obvi¬ 
ous  from  consideration  of  the  energy  balance  and  thermal  softening 
kinetics  in  (11):  larger  values  of  /?  lead  to  larger  amounts  of  energy 
converted  to  temperature  rise,  which  in  turn  leads  to  increased  ther¬ 
mal  softening.  Reasons  for  relatively  large  local  values  y?  —  1  in  the 
vicinity  of  regions  of  intense  plastic  strain  include  (i)  local  dislocation 
content  approaching  saturation  such  that  rates  of  hardening  and  en¬ 
ergy  storage  are  minimal;  (ii)  large  magnitudes  of  slip  rates  contrib¬ 
uting  to  the  denominator  in  (16);  and  (iii)  large  temperature  rise 
contributing  to  the  product  Odp/dO  <  0  in  (19).  Shear  localization 
due  to  adiabatic  heat  generation  and  thermal  softening  in  ductile 
metals  can  occur  similarly  in  compression,  shear,  and  mixed-mode 
loading  [36].  Shown  in  Fig.  1(e)  is  the  dislocation  density  associated 
with  cumulative  hardening  in  (13),  while  shown  in  Fig.  1(f)  is  the 
residual  elastic  volume  change  j  -  1  from  (15).  Dislocation  densities 
of  the  same  order  of  magnitude  have  been  predicted  elsewhere 
[24,25]  in  grain-scale  crystal  plasticity  simulations  of  FCC  metals. 
The  volume  changes  predicted  here  are  small  and  do  not  significantly 
affect  the  pressure  (e.g.,  Ap  a;  1.5  x  1CT4K  «  12  MPa),  and  occur  in 
regions  where  the  dislocation  density  is  highest.  According  to  (15), 
the  dislocation  density  and  the  residual  volume  change  are  linearly 
related,  with  (J  -  1)/pT  varying  with  temperature  as  a  result  of  the 
affect  of  temperature  on  the  elastic  moduli. 


In  summary,  the  results  in  Fig.  1  demonstrate  positive  correla¬ 
tions  among  temperature  rise,  plastic  strain,  and  the  heat  dissipa¬ 
tion  fraction  p.  Increases  in  parameter  Q,  via  processing  steps  to 
increase  the  capacity  for  stored  energy  of  mechanical  working, 
would  reduce  p,  heat  conversion,  and  the  tendency  for  localization 
of  plastic  deformation.  Texturing  could  conceivably  increase 
stress-carrying  capacity,  for  example,  via  inclusion  of  a  higher  pro¬ 
portion  of  the  stiffer  grains  evident  in  Fig.  1(a).  A  more  extensive 
study  would  consider  a  large  number  of  microstructure  morpholo¬ 
gies,  grain  orientation  sets,  and  loading  rates;  however,  the  present 
mesoscale  simulation  is  sufficient  to  demonstrate  mechanisms  to 
be  exploited  in  what  follows.  Specifically,  two  mechanisms  for  pos¬ 
sible  performance  enhancement  -  grain  orientation  and  energy 
storage  -  are  considered. 


4.  Coarse-grained  macroscopic  model 

Homogenized  results  of  a  number  of  polycrystal  calculations  are 
used  to  provide  parameters  for  a  macroscopic  plasticity  model 
used  in  impact  simulations.  Polycrystalline  aggregates  consisting 
of  300  crystals  are  subjected  to  uniform  shearing  at  a  constant  rate 
via  F 3  =  1  +  yt,  where  the  shear  strain  is  y  =  yt.  The  shear  bound¬ 
ary  conditions  at  an  applied  rate  of  y  =  104/s  are  those  deemed 
most  crucial  in  a  ballistic  perforation  event  [27],  Each  grain  in 
the  aggregate  is  sheared  adiabatically  under  the  same  total  defor¬ 
mation  gradient  [30,4],  and  each  grain  has  the  same  volume.  In  the 
polycrystal  calculations,  two  initial  textures  are  considered:  ran¬ 
dom  and  cubic.  The  random  texture  is  representative  of  the  simu¬ 
lation  in  Section  3  and  results  in  Fig.  1 ,  while  the  cubic  texture  can 
be  obtained  in  a  ductile  plate  by  rolling  [34].  In  a  previous  work  [6], 
a  number  of  different  rolled  textures  [34]  were  considered,  with 
cubic  found  to  give  the  greatest  increase  in  the  shear  stress  compo¬ 
nent  thought  to  provide  the  primary  resistance  to  plug  formation 
[19,27].  For  the  random  texture,  each  grain  is  assigned  an  indepen¬ 
dent  set  of  three  Euler  angles  with  a  random  number  generator, 
providing  for  isotropic  average  behavior.  To  obtain  shear  strength 
properties  appropriate  for  a  plate  with  a  cubic  texture  subject  to 
perforation,  grains  are  all  aligned  with  one  cube  axis  parallel  to 
the  normal  direction  of  the  A1  plate,  that  is,  parallel  to  x3,  the  direc¬ 
tion  of  perforation.  The  other  two  cube  axes  are  distributed  ran¬ 
domly  among  the  grains  in  the  plane  of  the  plate,  the  x1  -  x2 
plane.  This  provides  for  transversely  isotropic  symmetry,  and  al¬ 
lows  for  an  average  resistance  of  the  plate  to  perforation  by  a  cylin¬ 
drical  projectile  that  would  arise  from  a  true  cubic  texture  with  the 
three  cube  axes  aligned  with  the  normal,  transverse,  and  rolling 
directions,  respectively  [6[.  For  each  initial  texture  (i.e.,  random 
or  cubic),  values  of  Q  =  1 , 4,  or  8  are  also  considered,  to  investi¬ 
gate  possible  strengthening  effects  from  energy  storage. 

The  volume-averaged  shear  stress  of  the  aggregate, 
t  =  V-1  f  andV,  is  shown  in  Fig.  2.  Each  grain  can  support  a  differ¬ 
ent  individual  shear  stress  as  a  result  of  its  particular  lattice  orien¬ 
tation.  Cubic  texturing  provides  a  ~  25%  increase  in  peak  shear 
strength  over  aggregates  with  random  texture.  Increases  in  energy 
storage  factor  Q  result  in  a  slight  increase  in  shear  strength,  for 
example  for  random  grains  at  y  =  0.4,  f  =  0.446  GPa  for  Q  =  1 
compared  to  t  =  0.471  GPa  for  0  =  8.  The  strength  difference 
arises  from  the  greater  heat  generated  for  Q  =  1,  leading  to  more 
thermal  softening  relative  to  the  prediction  for  0  =  8.  Since  all 
grains  are  assigned  the  same  deformation,  no  representation  of 
localization  is  included  in  the  results  of  Fig.  2  that  would  enhance 
the  strength  distinction  afforded  by  different  values  of  O. 

The  average  dissipation  fraction  ~p  =  V~'  f  pdV  is  shown  in 
Fig.  3.  Each  grain  can  provide  a  different  value  of  the  local  dissipa¬ 
tion  fraction  (i.e.,  the  integrand)  as  a  result  of  its  own  orientation. 
In  all  cases,  this  quantity  decreases  initially  in  conjunction  with 
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Fig.  2.  Average  shear  stress  versus  shear  strain  for  A1  2139  polycrystals  of  various 
microstructures. 


work  hardening  and  dislocation  multiplication,  then  increases 
asymptotically  towards  unity  at  large  strains  as  the  dislocation 
content  saturates  [36,4].  Small  differences  in  the  dissipation  frac¬ 
tion  arise  between  random  and  cubic  textures.  Large  differences 
arise  from  the  choice  of  Q,  as  expected  from  (19).  As  noted  previ¬ 
ously,  commercially  pure  polycrystalline  A1  exhibits  an  average 
cumulative  value  of  fi  ~  0.93  [11].  Slightly  higher  values,  e.g. 
~  0.95  at  strain  levels  approaching  50%,  have  been  noted  for  pure 
aluminum  single  crystals  [10].  The  energy  storage  is  presumably 
greater  in  polycrystals  because  of  concentrated  residual  elastic  en¬ 
ergy  near  grain  boundaries  and  triple  points,  as  has  been  predicted 
in  mesoscale  simulations  of  FCC  polycrystals  [8].  Also  shown  in 
Fig.  3  are  high-rate  experimental  results  from  Rosakis  et  al.  [26] 
of  rate-insensitive  alloy  Al  2024-T351,  where  the  compressive 
plastic  strains  eP  reported  in  that  work  are  converted  to  equivalent 
shear  strains  using  y  a;  ^3eP.  The  adjustable  parameter  Q  dictates 
the  ratio  of  the  rate  stored  residual  elastic  energy  to  inelastic  stress 
power,  as  is  clear  from  Fig.  3.  The  larger  the  value  of  Q,  the  greater 
the  stored  energy.  Experimental  data  on  ~p  is  not  available  for  Al 
2139  which  is  also  fairly  rate  insensitive  ([6]  and  references  there¬ 
in),  but  the  results  in  Fig.  3  demonstrate  that  a  physically  reason¬ 
able  range  of  /)  is  provided  by  1  ^  Q  C  8.  The  dissipation  fraction 
of  Al  2024-T351  would  seem  to  correspond  to  Q  between  4  and 
8,  while  pure  Al  would  correspond  to  12  a;  1 .  On  the  other  hand, 
for  alloys  of  different  composition  or  processing  routes,  12  may 
deviate  from  these  values.  Differences  in  12  may  arise  among  differ¬ 
ent  alloys  because  of  variations  in  composition  and  microstructure. 
Quantitative  expressions  relating  12  to  precipitate  sizes  and  volume 
fractions,  impurity  concentrations,  and  dislocation  distributions 
remain  to  be  determined  by  careful  microscopic  experiments  and 
numerical  simulations  at  the  scale  of  discrete  dislocation  dynamics 
or  molecular  dynamics. 
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Fig.  3.  Average  heat  dissipation  fraction  versus  shear  strain  for  Al  2139  polycrystals 
of  various  microstructures;  experimental  data  points  estimated  from  Rosakis  et  al. 
[26]. 


The  plasticity  component  of  the  macroscopic  model  used  here 
to  represent  Al  2139  in  ballistic  simulations  follows  from  Zerilli 
and  Armstrong  [37].  For  FCC  metals  such  as  Al,  the  yield  surface 
in  terms  of  deviatoric  strength  is 

a  =  C0  +  C2£p  exp(-C30  +  C4dloge),  (21) 

where  s  =  ((2/3)D  :  D)1/2  is  the  effective  total  strain  rate  and  D  is 
the  symmetric  part  of  the  total  velocity  gradient  L.  Symbols 
Co,C2,C3,C4,  and  N  denote  fitting  parameters,  and  G  is  the  macro¬ 
scopic  shear  modulus.  Values  corresponding  to  the  Al  alloy  of  pres¬ 
ent  interest  with  different  textures  are  listed  in  Table  2;  specifically, 
C4  =  0  implies  rate  insensitivity,  and  C2  is  higher  for  the  cubic  tex¬ 
ture  than  the  random  texture,  leading  to  the  increased  flow  stress  of 
the  former  evident  in  Fig.  2.  These  parameters  are  obtained  by  fit¬ 
ting  (21)  to  the  results  of  the  polycrystal  simulations  of  Fig.  2,  again 
assuming  adiabatic  conditions.  Note  that  the  material  is  treated  in 
the  coarse  scale  model  as  isotropic,  a  rigorous  assertion  only  for 
Al  plates  of  the  random  texture.  The  cubic  orientation  would  afford 
transverse  isotropy,  as  discussed  above.  Thus,  the  normal  strength 
of  cubic-textured  plate  in  the  x3-direction  and  the  shear  stresses 
au  and  er23  may  be  over-  or  underestimated  by  the  assumption  of 
isotropy.  However,  the  predictions  of  dynamic  perforation  will  be 
accurate  so  long  as  the  primary  resistance  to  plugging  is  the  shear 
stress  t  =  cr13  tangential  to  the  penetration  direction  [19,27].  Also 
shown  in  Table  2  are  properties  for  iron  used  to  model  the  projectile 
in  impact  simulations  that  are  described  in  Section  5.  The  dynamic 
yield  stress  of  iron  is  given  by  the  BCC  Zerilli-Armstrong  model 
[37]: 

a  =  C0  +  C,  exp(-C30  +  C4dloge)  +  C5eP,  (22) 

with  Ci  and  C5  additional  parameters.  The  local  macroscopic  energy 
balance,  i.e.  the  analog  of  (8)  under  adiabatic  conditions,  is 

pci)  =  Jin  :  Dp  +  pp,0J),  (23) 

where  c  accounts  for  effects  of  the  equation  of  state  of  the  material 
that  relates  mass  density,  temperature,  and  defect  content  to  hydro¬ 
static  pressure.  The  effect  of  residual  elastic  volume  changes  in  (15) 
is  captured  in  a  modified  Mie-Gruniesen  equation  of  state  as 

p  =  (Kfl  +  K2p2  +  K3p?){  1  -  Tp/2)  +  rpocA0(l  +  p),  (24) 

where  K,l<2,  and  K3  are  dilatational  elastic  constants,  T  is  Grunie- 
sen’s  parameter,  and  p  =jj~'  -  1  measures  recoverable  elastic  vol¬ 
ume  changes.  The  reference  mass  density  is  p0  =  pj  =  p  detF. 
Library  property  values  for  aluminum  [18]  are  used  in  (24)  in  the 
ballistic  simulations  that  follow.  In  the  macroscopic  model  the 
residual  volume  change  is  found  by 

J=1  +  (Es/G)(dG/dp-G/K) 

=  1  +G-\dG/dp  -  G/K )  Dpdt.  (25) 

The  integral  in  the  final  term  of  (25)  is  the  cumulative  macroscopic 
stored  energy  per  unit  reference  volume  Es,  and  G  is  the  polycrystalline 


Table  2 

Parameters  for  coarse-scale  plasticity  model. 


Parameter 

Al  2139,  random 

Al  2139,  cubic 

Iron1 

C0  [GPa] 

0.40 

0.40 

0.65 

C,  [GPa] 

- 

- 

1.03 

C2  [GPa] 

2.5 

4.5 

- 

C3  [1/K] 

0.0045 

0.0045 

0.0070 

C4  [1/K] 

0 

0 

4.16(10)  4 

C5  [GPa] 

- 

- 

0.27 

N 

0.30 

0.30 

0.29 

G  [GPa] 

26.0 

26.0 

80.2 

1  Zerilli  and  Armstrong  [37], 
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shear  modulus  listed  in  Table  2.  The  first  relationship  in  (25)  is  recog¬ 
nizable  as  the  macroscopic  (polycrystalline)  analog  of  the  single  crystal 
expression  for  residual  volume  change  in  the  first  of  (15).  The  dissipa¬ 
tion  fraction  ji  in  (23)  and  (25)  is  treated  as  a  polynomial  function  of  the 
effective  plastic  strain: 

/f  -  -1  1  B  i£p  +  B2ej  +  B3Ep.  (26) 

Values  of  B i,  B2,  and  B3  are  calibrated  to  the  results  of  polycrystal 
simulations  of  Fig.  3  and  are  listed  in  Table  3  for  each  microstruc¬ 
ture.  The  projectile  is  assigned  the  same  properties  in  all  simula¬ 
tions,  with  the  default  value  of  /)  =  1  and  library  equation-of-state 
parameters  [15], 

5.  Terminal  ballistics  predictions 

Twelve  impact  simulations  are  performed  to  demonstrate  pos¬ 
sible  effects  of  initial  texture  and  energy  storage  on  perforation 
resistance  of  an  Al-Cu-Mg-Ag  plate.  These  include  plates  of  one 
of  six  different  A1  alloy  microstructures,  each  subjected  to  normal 
(i.e.,  null  obliquity)  penetration  by  an  iron  cylindrical  projectile 
at  an  impact  velocity  V\  of  either  1000  m/s  or  1500  m/s.  The 
dimensions  of  the  projectile  are  diameter  20  mm  and  length 
22  mm.  The  plate  is  25.4  mm  thick  and  254  mm  x  254  mm  square 
in  the  lateral  directions.  These  dimensions  are  sufficiently  large  to 
enable  the  projectile  to  penetrate  completely  through  the  target 
before  interactions  of  the  projectile  with  release  waves  from  the 
lateral  edges  occur.  Trends  in  results  reported  below  are  insensi¬ 
tive  to  refinement  of  the  mesh  via  a  ~  50%  increase  in  the  number 
of  elements. 

The  macroscopic  models  outlined  in  (21 )— (26),  with  properties 
listed  in  Tables  2  and  3,  are  used  for  the  plate  and  projectile.  Prop¬ 
erties  for  the  A1  plate  are  varied  according  to  microstructure  as 
listed  in  Tables  2  and  3:  two  textures  each  with  three  values  of  Q 
provide  a  total  of  six  property  sets  corresponding  to  six  different 
microstructures.  The  EPIC  dynamic  finite  element  code  [15]  is  used, 
in  which  elements  are  converted  to  interacting  particles  at  large 
effective  strains  (>  0.5)  to  avoid  numerical  difficulties  associated 
with  severe  mesh  distortion  [16,7].  In  all  simulations,  the  ductility 
of  the  aluminum  plate  is  assigned  a  value  of  unity;  when  strained 
beyond  this  point,  the  deviatoric  strength  of  the  material  degrades 
to  zero  locally.  Library  properties  [15]  for  ductility  and  failure  of 
the  iron  comprising  the  projectile  are  used  in  all  12  simulations. 

Representative  results  in  Fig.  4  depict  the  simulation  after  im¬ 
pact  at  1500  m/s,  with  properties  of  the  plate  corresponding  to  ran¬ 
dom  grain  orientations  and  Q  =  1.  The  plug  that  typically  emerges 
in  perforation  of  ductile  plates  impacted  at  high  velocity  [19,27]  is 
evident  in  Fig.  4.  Stress  concentrations  along  the  circumference  of 
the  perforated  zone  and  symmetrically  located  in  two  other  inte¬ 
rior  regions  of  the  plate  are  evident,  with  magnitudes  on  the  order 
of  the  maximum  yield  stress  of  the  material  (Fig.  2). 

Shown  in  Fig.  5  is  the  ratio  of  the  residual  velocity  VR,  i.e.  the 
average  steady-state  velocity  of  the  projectile  (in  the  direction  nor¬ 
mal  to  the  plate)  after  complete  perforation,  to  the  impact  velocity 
V/  for  all  12  simulations.  The  lower  the  value  of  VR/Vi,  the  more 
resistant  the  plate  to  dynamic  perforation.  For  each  microstructure, 
this  ratio  decreases  with  decreasing  impact  velocity.  Use  of  a  cubic 


Fig.  4.  Simulation  of  perforation  of  A1  2139  plate  (random  texture.fl  =  1)  by  iron 
cylinder  at  initial  velocity  of  1500  m/s,  shown  here  at  100  ps  after  initial  impact. 
Effective  stress  contours  in  aluminum  plate  also  shown. 


Fig.  5.  Normalized  residual  velocity  of  projectile  after  perforation  of  A1  plates  with 
inelastic  properties  corresponding  to  various  microstructures. 

texture  over  a  random  texture  seems  to  offer  a  ~10%  reduction  in 
Vr/Vi  at  an  impact  velocity  of  1 500  m/s  and  a  30%  or  greater  reduc¬ 
tion  at  an  impact  velocity  of  1000  m/s.  Recall  however  that  the 
assumption  of  isotropy  introduces  some  uncertainty  into  these 
predictions;  a  more  thorough  treatment  would  fit  the  polycrystal 
model  results  to  an  anisotropic  yield  surface  [20],  Residual  elastic 
energy  storage  in  the  plate  also  influences  performance,  with 
Vr/V(  reducing  by  ~5%  as  Q  is  increased  from  1  to  8  for  an  impact 
velocity  of  1500  m/s,  and  VR/V,  reducing  by  about  ~20%  as  Q  is  in¬ 
creased  from  1  to  8  in  the  random  microstructure  at  an  impact 
velocity  of  1000  m/s.  The  effect  of  stored  elastic  energy  is  even 
stronger  for  the  cubic  microstructure  at  an  impact  velocity  of 
1000  m/s,  wherein  VR/Vi  is  reduced  by  over  50%,  from  0.225  to 
0.105,  as  Q  is  increased  from  1  to  8.  The  improvement  in  perfor¬ 
mance  with  increasing  Q  occurs  because  as  less  adiabatic  heating 
takes  place  in  the  penetration  zone,  less  localized  deformation  re¬ 
sults,  delaying  failure  of  the  material. 


Table  3 

Parameters  for  dissipation  fraction  in  coarse-scale  model. 


Parameter 

Random  texture 

Cubic  texture 

0  =  1 

0  =  4 

0  =  8 

0  =  1 

0  =  4 

0  =  8 

B, 

-0.618 

-2.457 

-4.873 

-0.673 

-2.661 

-5.224 

b2 

0.861 

3.478 

7.054 

1.025 

4.101 

8.176 

b3 

-0.331 

-1.350 

-2.776 

-0.415 

-1.673 

-3.378 

6.  Conclusions 

A  crystal  plasticity  model  is  developed  for  precipitate-hardened 
aluminum  alloys  accounting  for  large  deformations,  thermody¬ 
namics,  nonlinear  elasticity,  and  energy  storage  mechanisms  and 
volume  changes  associated  with  residual  stress  fields  of  lattice  de¬ 
fects.  Direct  numerical  simulation  of  a  polycrystalline  microstruc¬ 
ture  with  polygonal  grains  demonstrates  variations  from  grain  to 
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grain  in  stress  fields  due  to  anisotropy,  as  well  as  positive  correla¬ 
tion  of  the  heat  dissipation  fraction,  adiabatic  temperature  rise, 
and  plastic  strain  localization.  A  macroscopic  plasticity  model  is 
calibrated  to  results  of  polycrystal  simulations.  The  macroscopic 
model  accounts,  via  choice  of  strain  hardening  parameters,  for  dif¬ 
ferences  in  shear  strength  afforded  by  texturing.  The  macroscopic 
model  also  captures  effects  of  stored  energy  that  depends  on  the 
strain  history,  and  the  model  accounts  for  volume  changes  from 
defect  generation  in  the  equation  of  state  for  the  hydrostatic  pres¬ 
sure.  Impact  simulations  demonstrate  that  improved  resistance  to 
ballistic  perforation  may  be  afforded  by  texturing  to  increase  shear 
strength  or  by  processing  steps  that  increase  the  activity  of  dy¬ 
namic  energy  storage  mechanisms  relative  to  dissipation  mecha¬ 
nisms.  In  the  future,  in  addition  to  shear  localization,  fracture 
should  be  given  more  consideration.  For  example,  mesoscale  sim¬ 
ulations  of  the  mechanism  of  spall  fracture  [33]  that  may  occur 
in  blast  loading  may  provide  insight  into  effects  of  statistical  vari¬ 
ability  of  microstructures,  e.g.  heterogeneous  spatial  distributions 
of  grains  and  inclusions,  on  resistance  to  failure. 
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B  FORCH 
SKARNA 
J  MCCAULEY 
RDRL  WMS 
T  JONES 
RDRL  WMM  A 
J  ANDZELM 
W  NOTHWANG 
RDRL  WMM  D 
B  CHEESEMAN 
G  GAZONAS 
RDRL  WMT  A 
R  DONEY 
S  SCHOENFELD 
RDRL  WMT  C 
T  BJERKE 
T  FARRAND 
M  FERMEN-COKER 
S  SEGLETES 
B  SCHUSTER 
A  TANK 
W  WALTERS 
C  WILLIAMS 
RDR  WMT  D 
S  BILYK 
D  CASEM 

J  CLAYTON  (10  CPS) 
D  DANDEKAR 
N  GNIAZDOWSKI 
M  GREENFIELD 
Y  HUANG 
R  KRAFT 
B  LOVE 


